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© The invention relates to an apparatus and meth- 
ods for generating multiparameter optical data that 
characterize a population of cells. The invention in- 
cludes the steps of scanning the cell population with 
a beam to produce sets of digital data samples, 
each sample set of digital data representing mul- 
tiparameter optical interactions from a specific loca- 
tion within the cell population; storing the digital 
data, e.g., in a computer memory; locating a cell 
within the population, e.g., by comparing the digital 
data to a preselected threshold value; defining a 
neighborhood around the digital data representing 
the located cell; estimating a background level for 
the neighborhood based upon digital data corre- 
sponding to locations outside the neighborhood; and 
correcting each of the samples corresponding to the 
neighborhood with the estimated neighborhood 
background level to generate the optical data. The 
invention further relates to specific methods of back- 
ground correction and data calibration as well as 
specific sampling features to enable precise esti- 
mates of multiple cellular constituents or other cell 
properties at high rates of speed. 
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METHOD AND APPARATUS FOR MEASURING OPTICAL PROPERTIES OF BIOLOGICAL SPECIMENS 



This invention relates to measuring multiple 
optical properties of biological specimens, such as 
a population of cells, e.g., blood cells, at high rates 
of speed using computer controlled instruments 
that scan the biological specimens. 

Examples of such computer controlled instru- 
ments include flow cytometers, automated blood 
cell analyzers and blood cell differential classifiers. 

A flow cytometer is an instrument that 
hydrodynamicaliy focuses a fluid suspension of 
cells into a single file stream that passes through 
an examination zone. A focused light beam illumi- 
nates the cells in this zone and the instrument 
measures optical interactions of the light with the 
cells such as, for example, multiple wavelength 
absorption, scatter as a function of angle, and 
fluorescence as a function of either wavelength or 
polarization. This type of instrument permits the 
study of living cells in addition to those which have 
been chemically treated, for example, by staining. 
Flow cytometry techniques certain constituents or 
structures.particularly those present on the cell sur- 
face, to be quantitatively characterized at cell rates 
of a thousand or more cells per second. 

Blood cell analyzers typically consist of a com- 
puterized microscope that automatically classifies 
various types of white blood cells and flags and 
counts all abnormal cells in a specimen. Using 
such an automated instrument, the operator can 
also view the cells manually, stop the analyzer 
temporarily for making visual morphological ob- 
servations, or review abnormal cells in greater de- 
tail. During the automatic count, the X-Y coordi- 
nates of every encountered abnormal blood cell are 
stored. Therefore, the operator has the option of 
using a review mode in which abnormal cells are 
automatically and individually acquired and fo- 
cused, and can view these cells on a T.V. monitor 
or through a binocular microscope. 

A blood cell differential classifier typically con- 
sists of a computer-controlled microscope having a 
stage that is driven by stepper motors. A light 
source, such as a xenon arc lamp, illuminates cells 
and the classifier uses various sensors such as a 
silicon photodiode array to measure the optical 
interaction between the cells and the light. All of 
the cells in a given area are illuminated at the 
same time. 

We describe below a method and apparatus for 
generating optical data that accurately estimates 
multiple constituents and simultaneously character- 
izes a number of morphological properties of each 
of a population of cells. This method includes the 
steps of, and means for scanning the cell popula- 
tion with a beam to produce digital data samples, 



the different digital data samples representing mul- 
tiple optical measurements at different locations 
within the cell population; storing the digital data, 
e.g., in a computer memory; locating a cell within 

5 the population, for example by comparing digital 
data derived from the stored digital data to a 
preselected threshold value; defining a neighbor- 
hood around the located cell; estimating a back- 
ground level or individual background levels for all 

w sample points in the neighborhood based upon 
stored digital data corresponding to locations out- 
side the neighborhood; and correcting each of the 
digital data samples corresponding to the neighbor- 
hood with the estimated neighborhood background 

15 level to generate the optical data. The beam used 
in the invention is of electromagnetic radiation, e.g., 
laser light, X-rays, or infrared radiation. 

The sample correcting step may be performed 
by subtracting the corresponding background level 

20 from each of the stored digital data samples in a 
neighborhood. 

The term "digital data sample" describes digi- 
tal information that may be stored in a computer 
memory and that is sampled from an analog optical 

25 signal from the population of cells in a specimen 
when scanned with the beam. The optical data 
results from correction of the digital data for back- 
ground. Further processing of the optical data gen- 
erates corresponding optical and morphological 

30 property values of the population of cells that may 
be used for precise measurements of cellular con- 
- stituents and as representations of cellular mor- 
phology. 

The term "constituent" as used herein refers to 

35 a specific part or component of a cell such as, for 
example, the DNA or HHA. The invention is ca- 
pable of determining the amount of a specific con- 
stituent in a cell or cell population with an accuracy 
of within about a few percent of the actual amount. 

40 The term "neighborhood" means a specific 

region of sample points that is created by the 
operator to encompass the digital data and that is 
set to statistically correspond to a single cell. That 
is, on average, only one cell will be located within a 

45 neighborhood. The size of such a neighborhood 
may be varied by the operator depending upon the 
size and concentration of the cells in the specimen. 
The neighborhood is used to allow an accurate 
estimate of cell constituents by determining the 

so background within each neighborhood, i.e., individ- 
ually for each cell. This allows for greater accuracy 
than a background level used for a number of 
different cells, because the background level varies 
across the examining surface, such as a micro- 
scope slide. 
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We also describe a method and apparatus for 
generating optical data that characterizes a 
population of cells that includes the steps of, and 
means for, scanning the cell population with a 
beam with a spot size comparable in size to the 
cells to be scanned to produce digital data sam- 
ples, in which different samples represent different 
locations within the cell population and the rate of 
sampling is such that the distance traveled per 
time as the spot is scanned between consecutive 
sampled locations approximates the spot size; and 
then processing the digital data to generate optical 
data that characterizes the cells in the population. 
The term "comparable" as used herein means 
within one order of magnitude of the size of the 
cells to be scanned. 

In general, the number of sample points re- 
quired to represent a cell is equal to the number of 
samples taken along a scan times the number of 
scans. When a small spot size of, e.g., one-half 
micron is used, as in the prior art, on the order of 
four hundred samples are required to accurately 
represent the oonstituents in a 10 micron by 10 
micron cell area. We can now have a spot size 
which is on the order of 10 microns in diameter 
and is sampled at a rate such that the distance 
traveled per time as the spot is scanned between 
consecutive sampled locations is approximately the 
spot size. This gives the benefit of processing a 
sample about four hundred times faster than with 
the small spot size without losing accuracy. As the 
size is decreased and the sampling rate is in- 
creased this benefit correspondingly decreases. 

This method may also include the steps of 
multiplying the neighborhood value by a velocity 
normalization factor proportional to the velocity at 
which the beam is scanned past each location to 
generate optical data calibrated for scanning ve- 
locity variations. 

We also describe a method and apparatus for 
generating calibrated optical data that characterizes 
a population of cells including the steps of, and 
means for, measuring the resulting intensity of a 
beam having a predetermined incident illumination 
intensity used to scan the cell population as a 
function of location of each sample within the scan; 
scanning the cell population with the beam to pro- 
duce digital data samples, different samples repre- 
senting different locations within the cell population; 
storing the digital data, e.g., in memory; normaliz- 
ing the stored digital data using the beam intensity 
measurements of the location of each sample with- 
in the scan; and then processing the normalized 
digital data to generate calibrated optical data that 
characterizes the cells in the population. The re- 
sulting beam intensity may be measured, for exam- 
ple, with a sensor, such as a photomultiplier. 

We also describe a method and apparatus for 



generating calibrated optical data that characterizes 
a population of cells, that includes the steps of and 
means for, scanning the cell population on a sur- 
face with a beam to produce digital data samples, 

5 the beam having a given illumination intensity, and 
different samples representing different locations 
within the cell population; storing the digital data, 
e.g., in a memory; scanning a region of the surface 
without cells to generate intensity calibration data; 

w generating an intensity normalization factor from 
the intensity calibration data and storing the inten- 
sity normalization factor in memory; locating a cell 
within the population, e.g., by comparing the digital 
data to a preselected threshold; defining a neigh- 

75 borhood around the located cell, the neighborhood 
containing optical data derived from the samples of 
digital data within the neighborhood; and correcting 
the optical data in a neighborhood, e.g., by mul- 
tiplying the optical data by the intensity normaliza- 

20 tion factor, to generate the calibrated optical data. 

The surface without cells may be a calibration 
slide including a uniform fluorescent dye surface, a 
dye-filled cuvette or a uniform scattering surface. 
In another embodiment, the processing step 

25 may include multiplying the digital data by a ve- 
locity normalization factor proportional to the ve- 
locity at which the light beam is scanned by each 
location to generate optical data calibrated for 
scanning velocity variations. 

30 The method may also include the further steps 

of estimating background levels for the individual 
neighborhoods based upon stored digial data cor- 
responding to each optical parameter at locations 
outside that neighborhood; and correcting each of 

35 the digital data samples corresponding to the 
neighborhood to generate optical data, e.g., by 
subtracting the corresponding neighborhood back- 
ground level from each of the samples in the 
neighborhood. 

40 For each of the above methods, the neighbor- 

hood background levels can be determined by 
finding minimum digital data sample values in 
memory locations corresponding to actual locations 
adjacent the neighborhood for each optical param- 

45 eter of each sample in the neighborhood, or by 
averaging a plurality of digital data sample values 
in memory outside the neighborhood. Typically, the 
term "adjacent" means within five scans from the 
bounds of a neighborhood. 

so in a further embodiment, time corresponding 

approximately to the time the digital data was ac- 
quired is recorded and stored as synchronous digi- 
tal time data, whereby each sample has a cor- 
responding digital time point. As used herein, the 

55 terms "time corresponding approximately to the 
time the digital data was acquired" mean the time 
a given strip of scans is completed. The cell popu- 
lation may also be located on a movable surface, 
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the scanning may proceed in synchrony with a 
stepping motion of the surface, and the location 
corresponding to a sample may be derived from 
digital location data that may be generated, where- 
by each sample has a corresponding digital time 
point and location. The movable surface may be a 
microscope slide, a cuvette into which the cells are 
inserted, or any other surface upon which or into 
which cells may be deposited. The cells on the 
movable surface of all of the above methods and 
apparatus can be visually observed. 

In another embodiment, an additional step of 
rescanning the population of cells, generating new 
digital data, location data and time data, and com- 
paring the new data with the corresponding data 
from a previous scan to determine changes in the 
cell population may be added. 

In a further embodiment of the above methods 
and apparatus, a cell can be located by detecting a 
digital data sample corresponding to the edge of 
the cell first contacted by the beam. In another 
embodiment, the neighborhood is initially centered 
around the first edge sample. The method may 
further include the steps of determining the maxi- 
mum digital data sample within the neighborhood 
and then recentering the neighborhood with respect 
to the maximum sample before estimating a back- 
ground level for the neighborhood. 

We also describe a method and apparatus for 
measuring multiple cellular properties including the 
steps of, and means for, scanning the cell popula- 
tion with a beam to produce sets of digital data 
samples, different samples representing different 
locations within the cell population, e.g., storing the 
digital data in a memory; locating a cell within the 
population, e.g., by comparing the digital data to a 
set of preselected thresholds; defining a neighbor- 
hood around the located cell; summing all samples 
of optical data derived from the sets of digital data 
in a neighborhood to generate optical neighbor- 
hood values for the parameters for that neighbor- 
hood; comparing each parameter sample in a 
neighborhood with a set of preselected thresholds 
to generate a binary neighborhood pattern of sam- 
ples either above or below th thresholds for each 
parameter; storing the optical neighborhood values 
and binary neighborhood patterns in memory; and 
generating optical properties corresponding to the 
optical neighborhood values and morphological 
properties corresponding to the binary neighbor- 
hood patterns for each cell. 

The optical properties may be generated by 
applying math functions to one or more of the 
optical neighborhood values to generate optical 
property values of the cell corresponding to the 
neighborhood. The morphological properties may 
be generated by applying boolean functions to one 
or more of the binary neighborhood patterns to 



generate morphological property values of the cell 
corresponding to the neighborhood. 

The term "optical neighborhood value" de- 
scribes the sum of all optical data in a given 

5 neighborhood for one parameter. This optical 
neighborhood value may then be used as is or 
processed further to generate optical property val- 
ues proportional to specific cell constituents cor- 
responding to cells of the population. The term 

70 "binary neighborhood pattern" describes a word 
pattern in memory corresponding to a single neigh- 
borhood, i.e., cell, and a specific parameter, that is 
generated by comparing each sample point in a 
neighborhood to a threshold and registering a "1" 

75 if the point is above a given threshold or a "0" if 
not. By testing each point, a memory word of ones 
and zeroes is generated which is the binary neigh- 
borhood pattern. This pattern can be further pro- 
cessed to generate morphological property values 

20 of the cell. 

Each of these methods can include the further 
steps of estimating a set of background levels for 
the neighborhood based upon digital data corre- 
sponding to locations outside that neighborhood; 

25 and subtracting the neighborhood background lev- 
els from each of the samples of digital data cor- 
responding to the neighborhood to generate optical 
data. 

We also describe an apparatus for measuring 

30 multiple cellular properties that includes a beam 
source for generating a beam of electromagnetic 
radiation, e.g., light, X-rays or infrared radiation that 
may be controlled by the computer during the 
scanning process to illuminate the cells with a 

35 beam spot comparable in size to the cells to be 
scanned, an optical signal being generated as a 
result of the spot illuminating the cells; a surface 
upon which the cells are located; an optical path for 
directing the spot from the beam source to the 

40 cells on the surface; a scanner that is interposed in 
the optical path between the light source and the 
cells for scanning the spot across the surface so 
that the spot passes over a plurality of cells on the 
surface; one or more sensors for measuring the 

45 optical signals; an analog-to-digital converter ar- 
ranged to sample the optical signals at a specified 
rate and to produce digital data; whereby the sam- 
pling rate is set such that the distance traveled per 
time as the spot is scanned between consecutive 

50 sampled locations approximates the spot size; data 
processors for correcting the digital data in the 
memory, generating optical property values from 
corrected data points corresponding to multiple 
constituents of individual cells, and processing the 

55 values to determine multiple cellular morphological 
properties; and a memory for storing the digital 
data produced by the converters. 

In further embodiments, the spot size used in 
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the apparatus and in all of the methods according 
to the invention is 3 to 50 microns or may be 
adjusted to be comparable in size to the size of the 
cells to be scanned. In addition, the wavelength 
distribution of the beam may be controlled by the 
data processor. 

The scanner of the apparatus 
may be a resonant galvanometer scanner. Further- 
more, the surface upon which the ceils may be 
located may be a microscope slide, a cuvette into 
which the cells are inserted, or any other surface 
upon which or into which the cells may be depos- 
ited. 

In another embodiment of the above apparatus 
and methods, the cell population may be on a 
movable surface, the scanning may proceed in 
synchrony with a stepping motion of the surface 
and digital location data may be generated from 
which the location corresponding to a sample can 
be derived. 

The digital data produced by the cells and the 
beam in any of the above methods and apparatus 
may be derived from sets of analog optical signals. 
Furthermore, these optical signals may be fluores- 
cence, forward angle light scatter, extinction, or 
wide angle light scatter. 

In another arrangement, we describe a method 
for generating calibrated optical data that character- 
izes a population of cells in which a number of 
cells of that population has a constant optical value 
for a given parameter. This method corrects for 
beam illumination and light collection differences 
along each scan line. This method includes the 
steps of scanning the cell population with a beam 
along a scan line to produce digital data samples, 
different samples representing different locations 
along the scan line within the cell population; stor- 
ing the digital data; locating a cell within this popu- 
lation; defining a neighborhood around the located 
cell, the neighborhood containing optical data de- 
rived from the stored digital data samples within 
the neighborhood; determining the location of the 
maximum digital data sample within the neighbor- 
hood and then recentering the neighborhood with 
respect to the maximum sample; summing all op- 
tical data values in a recentered neighborhood to 
generate an optical neighborhood value for that 
neighborhood; determining the most frequently oc- 
curring optical neighborhood value to define a 
mode optical value; selecting a subpopulation of 
cells with optical neighborhood values within a pre- 
determined range around the mode optical value; 
determining a set of average optical neighborhood 
values of cells in the subpopulation as a function of 
cell location along the scan line; computing an 
array of correction coefficients equal to the ratio of 
the mode optical value to the average optical value 
for each location along the scan line to generate an 



illumination normalization factor; and correcting op- 
tical data for each cell at different locations along 
the scan line by multiplying the optical neighbor- 
hood value for that cell by the corresponding illu- 

5 mination normalization factor. The term "mode op- 
tical value" describes the optical neighborhood val- 
ue that occurs most frequently in a complete sam- 
ple run. The term "illumination normalization fac- 
tor" describes an array of correction coefficients 

10 equal to the ratio of the mode optical value to the 
average optical neighborhood value of all cells for 
each location along the scan line. 

In another arrahgement, we describe a method 
for generating calibrated optical data that character- 

15 izes a population of cells. The method includes the 
steps of scanning a population of calibration par- 
ticles comprising a known constant parameter with 
a beam along a scan line to produce digital calibra- 
tion data samples, different calibration samples re- 

20 presenting different locations along the scan line; 
storing the digital calibration data; locating a cali- 
bration particle within the population; defining a 
neighborhood around he located particle, the 
neighborhood containing optical calibration data de- 

25 rived from the stored digital calibration data sam- 
ples within the neighborhood; determining the loca- 
tion of the maximum digital calibration data sample 
within the neighborhood and then recentering the 
neighborhood with respect to the maximum digital 

30 calibration data sample; summing all optical cali- 
bration data values in a recentered neighborhood to 
generate an optical neighborhood calibration value 
for that neighborhood; determining the most fre- 
quently occurring optical calibration neighborhood 

35 value to define a mode optical calibration value; 
determining a set of average optical calibration 
neighborhood values as a function of particle loca- 
tion along the scan line; computing an array of 
correction coefficients equal to the ratio of the 

40 mode optical calibration value to the average op- 
tical calibration value for each location along the 
scan line to generate an illumination normalization 
factor; scanning the cell population with a beam 
along the scan line to produce digital data samples, 

45 different samples representing different locations 
along the scan line within the cell population; stor- 
ing the digital data; deriving optical data from the 
stored digital data; and correcting optical data for 
each cell at different locations along the scan line 

so by multiplying the optical data for that cell by the 
corresponding illumination normalization factor. 

Other features and advantages of the invention 
will be apparent from the following description of 
the preferred embodiments, and from the claims. 

55 The drawings will first briefly be described. 

FIG. 1 is a schematic diagram of an instrument 
for measuring multiple optical properties of bio- 
logical specimens at high rates of speed. 
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FIG. 2 is a block diagram of the electromechan- 
ical circuit used in the instrument shown in Fig. 
1. 

FIG. 3 is a flow chart of the general optical 
signal data acquisition loop. 
FIG. 4 is a schematic drawing of the scan pat- 
tern of the laser beam. 

FIGS. 5a and 5b are graphs illustrating the prop- 
er sampling rate according to the invention. 
FIG. 6 is a series of graphs representing analog 
optical data and the corresponding digital data. 
FIG. 7 is a flow chart of the data processing 
function steps used to manipulate data stored in 
memory. 

FIG. 8 is a monitor screen display for an actual 
test run charting red fluorescence value (x-axis) 
against forward angle scatter value (y-axis). 
Fig. 9 is a flow chart of the data processing 
steps of a calibration method. 



Mechanical and Optical Systems 

Referring to Fig. 1, the instrument 10 includes 
a light source 12, a mirror scanner 14, such as a 
resonant galvanometer scanner, an epi-illumination 
microscope 16, a stepper motor controlled stage 
18, light detectors 20, 22, 24, 26, and various 
associated optical components which will be de- 
scribed below. Light source 12 produces a light 
beam 12a that reflects off of scanner 14 to produce 
a scan beam 12b and finally illuminates a scan 
spot i2c of a fixed diameter or size on a specimen 
plane or surface 28. Specimen plane or surface 28 
is positioned on stage 18. Light source 12 is a 
laser such as, for example, a Helium Neon, Argon 
ion, Helium Cadmium, or solid state laser, depend- 
ing on the application. More than one laser may be 
used for a given application. If this is the case, the 
beams can be combined using a dichroic mirror so 
that they are coaxial. For some applications it may 
be desirable to control the intensity of the laser 
beams or shutter them under control of a com- 
puter. Lasers with multiple wavelength outputs may 
also be used, in which case, it may be desirable to 
use a computer-controlled filter, prism or Bragg cell 
to select a specific wavelength. 

After passing through a dichroic mirror 30, the 
laser beam is imaged by two lenses 32 and 34 
onto an epi-illumination* field stop 36 of microscope 
16. Resonant scanner 14 is located between lenses 
32 and 34 and scans the beam across the field 
stop when electrically driven. The focal lengths of 
lenses 32 and 34 and the deflection angle of scan- 
ner 14, which is proportional to the galvanometer 
drive voltage, control the size of the spot and the 
length of the scan at field stop 36 and thus, at 
specimen surface 28. The scanner is driven at 800 



Hz and the spot size at the specimen is 10 ix and 
the scan length at the specimen is 120 u. These 
are nominal values and can be changed by the 
user by rotating microscope nosepiece 38 bearing 

5 0b]ectives 39 from the nominal 20X to other higher 
magnifications to reduce the spot size and scan 
length or lower magnifications to increase them. 

Epi-illumination is used to illuminate the speci- 
men and to transmit fluorescent or scattered light 

ro to the viewing eyepiece 40 or to a film or video 
camera 42. The light transferring assembly 44 may 
contain a dichroic or partially or fully silvered mirror 
as well as an optical filter in the viewing path. 
These assemblies are interchangeable and the mi- 

15 croscope used in the described embodiment in- 
cludes a movable rod to exchange these assem- 
blies. 

The specimen surface 28 may be a slide upon 
which a tissue or cytology specimen is mounted or 

20 a cuvette into which a cell specimen is injected. 
Furthermore, microscope 16 may be an inverted 
microscope and specimen surface 28 may be a 
container or dish containing living cells. The optical 
instrument is designed to accept either a standard 

25 or inverted microscope. 

The cells of the specimen may be stained to 
alter their light scattering characteristics or cause 
them to fluoresce. The scatter or fluorescence of 
each cell is related either to 1) the amounts of 

30 specific cellular constituents such as DNA or pro- 
teins or 2) the amounts of one or more fluorescinat- 
ed antibodies which bind only to specific cell anti- 
gens. 

The position of the specimen with respect to 

35 scanning beam 12b is controlled by X and Y step- 
ping motors 46, 47. These motors are driven in 
steps that are of a size ranging from a fraction of 
the spot size to the spot size. This step size is 
under the control of a computer program. The 

40 stage is also provided with sensors to provide 
signals to the computer to indicate a "Home" or 
reference position io for the stage and to limit its 
travel. By always moving the stage to Home at the 
beginning of each run or periodically during a run 

45 (i.e., a group of scan strips), it is possible to obtain 
the absolute position of given cells on the speci- 
men surface so they can be manually reviewed or 
remeasured. For normal use, the stage can be 
moved manually by knobs attached to the motors. 

so The stage may also be manually controlled by 
buttons, a joystick or a computer mouse with a 
specifically programmed display. The microscope 
16 is focused manually by knob 62. 

Light scattered from the specimen in the for- 

55 ward direction is primarily related to cell size but 
may be altered by staining the cell. This scattered 
light passes through the microscope condenser 
lens 48, through a dichroic or partially silvered 
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mirror 50, and is imaged onto two sensors 20 that 
are electrically connected but spaced to allow the 
incident scan beam image to pass between them. 
Therefore, only light scattered by the specimen in 
the forward direction away from the axial light path 
is collected by sensors 20. The microscope can be 
used normally to illuminate a specimen by the 
standard illuminator 52 because the partially sil- 
vered or dichroic mirror 50, if suitably chosen, will 
not significantly reduce illumination. Mirror 50 may 
be placed into an assembly that mounts directly 
onto the microscope substage illuminator and is 
smail enough so that it does not interfere with 
focusing the condenser. The signals from sensors 
20 are amplified and become one input for the data 
acquisition circuit described in more detail below. 
Although not shown, an additional sensor may be 
placed behind or below sensors 20 to collect light 
from the incident scan beam in order to measure 
extinction. 

Light scattered at 90 degrees from the cells is 
useful for distinguishing certain types of ceils such 
as blood granulocytes, because cells with granules 
scatter light at oblique angles. Obliquely scattered 
light will be trapped, with the specimen slide or 
cuvette acting as a light pipe, and can be detected 
by a sensor 22 placed at the slide or cuvette edge. 
This sensor is mounted on the stage in contact 
with the slide or cuvette and its amplified signal is 
used as a second input for the data acquisition 
circuit. Again, this sensor does not interfere with 
normal microscope usage. 

Backscattered or 180 degree scattered light 
and fluorescent light are collected by objective lens 
38 at high numerical aperture and imaged back 
through the lenses 35, 34, and 32 to dichroic mirror 
54 which is designed to reflect the longer 
wavelength fluorescence if backscatter is not to be 
measured or to reflect some of the laser source 
wavelength if backscatter is measured. Dichroic 
mirror 30 transmits almost all of the laser 
wavelen.gth. Dichroic mirror 54 splits the light into 
two parts so as to measure backscatter at the laser 
wavelength and longer wavelength fluorescence, or 
two different wavelengths of fluorescence. Mirror 
54 can be partially silvered or reflect and transmit 
light based on polarization if applications require 
measuring fluorescence depolarization. Mirror 54 
reflects part of the incident energy through filter 56 
of the appropriate bandpass wavelength onto 
photomultiplier 24. Similarly, mirror 58 reflects the 
remaining energy through filter 60 to photomul- 
tiplier 26. The signals from the two photomultipliers 
are amplified and become additional inputs for the 
data acquisition circuit. Apertures have not been 
shown in the figures and are used in various places 
to reduce stray light reflected from the surfaces of 
the lenses and mirrors. 



Electromechanical System 

Referring to Fig. 2, the electromechanical sys- 
tem provides means of inputting signals from the 

5 light sensors 20, 22, 24 and 26, the scanner driver 
70, and the microscope stage 18, to the computer 
81 and outputting signals from the computer to the 
X stepper motor 46 and Y stepper motor 47 which 
move the microscope stage 18. This is accom- 

w plished by a commercially available circuit board 
80 which accepts four analog voltages, digitizes 
them at rates up to 100,000 Hz with analog-to- 
digital converter 80a and causes the resultant com- 
puter words to be stored in the computer memory 

75 under direct memory access (dma) control. Board 
80 also accepts digital values from inputs 82-90, 
which provide limit information from the stage, and 
provides digital output values on lines 92-104, 
which control the stage. It also controls the values 

20 of two analog voltages used to control the supply 
voltages to photomultipliers (light detectors) 24 and 
26. A detailed description follows. 

In one embodiment, circuit board 80 is a Data 
Translation Model 2828 (Marlborough, New York) 

25 which plugs into one of the slots on an 80286 IBM 
compatible computer. The computer is represented 
in Fig. 2 by block 81 . 

Furthermore, photomultipliers 24 and 26 are 
used to detect fluorescence or backscatter and 

30 photosensors 20 and 22 are used to detect forward 
and oblique scatter emissions from the specimen 
and are each connected to respective amplifiers 
76, 72, 78 and 74. These amplifiers have appro- 
priate gains to provide signal levels of -10 to +10 

35 volts to circuit board 80 at analog signal inputs 106, 
108, 110 and 112. The circuits 72 to 76 are D.C. 
operational or instrument amplifiers with D.C. "no 
signal" levels near -10 volts. As will be described 
below, it is necessary to synchronize the position 

40 of the scanning mirror 14 with the data acquired by 
the sensors and converted to a stream of digital 
data. The digital data stream may be stored as two 
blocks of 64,000 words each. For synchronization 
purposes, the D.C. level of circuit 78 is initially set 

45 to 0 volts, so that a negative synchronization pulse 
is easily detected. 

Synchronization is accomplished by using a 
square wave pulse synchronization signal gener- 
ated by the scanner mirror driver 70 which controls 

so the motion of scan mirror 14 through scanner 72. 
This signal is modified to a pulse by circuit 114 
and, in one embodiment, is added to the sensor 
signal from sensor 20. Although this reduces the 
digitization resolution of this sensor input by two, a 

55 separate data channel for synchronization is not 
required. Of course, the synchronization signal may 
be used as a separate input signal. The signal at 
input 112 is the forward scatter signal of sensor 20 
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and the pulse signal is negatively added near one 
scan extreme (as shown in Fig. 6). As will be seen 
below, this negative pulse is detected by the pro- 
gram and used to properly synchronize the digital 
data stored in the computer memory. The sampling 
rate is set by the user through an initialization 
program which allows the user to define a Protocol 
for each test. The 

Protocol is a monitor screen that the operator 
uses to set the sampling rate and the various test 
parameters, area scanned, threshold settings, etc. 
The number of parameters digitized is also preset 
and the amplifier gain settings and input/output 
relationship, i.e., linear or logarithmic, may be used 
as additional parameters. 

The levels of the digital outputs 92, 94, 96, 98, 
100, 102 and 104, of circuit board 80 are under the 
control of the computer program. The digital inputs 
82, 84, 86, 88, and 90 are read at specific times 
also determined by the control of the program. 
These outputs and inputs are used to control the 
movement of the .microscope stage 18 via X step- 
per motor 46 and Y stepper motor 47, which are 
each driven by translator circuits 116 and 117, 
respectively. The microscope stage is provided 
with limit switches which indicate when the stage 
has reached its limit of travel in the x and y 
directions. These switches generate signals on 
lines 18a-18d which are used as inputs 82 to 88, 
respectively, of board 80. Input 90 senses the 
pulses from scanner 70. Although not shown here, 
additional digital outputs may be used to control 
the wavelength of the light source by controlling 
specific light sources, shutters, or filter positions. 

The program controlled stage motion is de- 
signed to perform the following sequence that is 
depicted in the flow chart of Fig. 3. First, when the 
user initiates a test, both stepper motors are 
moved to a specific "Home" position. This is ac- 
complished by calling a program subroutine to 
pulse lines 94 and 100 on and off until inputs 82 
and 88 indicate that the stage has reached Home. 
An OR rate 118 passes the signal directly from 
output 94 to input 116a. Under program control, 
outputs 96 and 1 02 are set to produce the proper 
stage direction by producing signals received at 
inputs 116b and 117b, As soon as the stage 
reaches Home, the Y stepper is pulsed to move the 
stage 18 to the initial y position, then a subroutine 
is called to move the stage to the right in the x 
direction by changing the signal on output 96. In 
one embodiment, the pulse rate on output 94 is 
ramped-up in rate from about 100 up to 800 pulses 
per second (pps) for a fixed total number of pulses 
or distance, typically 100 pulses or steps. This is 
the ramp-up numbers This fixed ramp-up may also 
be adjusted by the program or by a commercially 
available ramp-up controller circuit (Metrobyte, 



Taunton, Massachusetts). 

The program produces a rate of 800 pps, 
which is the scanner frequency, at the end of the 
ramp-up, so the stepper is moving at full stepping 

5 speed at this time. The step size may be either 5 
or 10 microns or some other size depending on the 
level at output 98 which controls the translator 116 
through input 116c to produce full or half steps of 
motor 46. Thereafter, output 92 is programmed to 

10 provide an input to AND gate 120 to allow the 
stepper motor to be directly driven by the synchro- 
nization output of mirror scanner driver 70 stepping 
stage 18 to the right in the x direction at 800 steps 
per second, e.g., one step for each scan cycle. The 

75 stepping motion in conjunction with the scanning 
motion generated by scanner driver 70 which is 
perpendicular to the stepping motion creates the 
scan pattern shown in Fig. 4. In Fig. 4, the scan 
starts at the left, at the Home position, and forward 

20 and reverse scans of length "L" are produced until 
the end of one scan strip is reached. Such a strip 
typically encompasses 2500 forward and 2500 re- 
verse scans. An additional circuit may be provided 
between driver 70 and translator 116 to allow other 

25 ratios of scanner rate to stepper rate under control 
of the user through the Protocol. 

An additional parameter of the protocol is the X 
scan distance. This distance determines the length 
of one scan strip. This length can be used to 

30 calculate the size required for the data buffers by 
multiplying the number of parameters measured by 
the total number of data values digitized per scan 
strip. The board 80 digitizes inputs, stores them in 
a buffer, and returns a flag when the buffer is full. 

35 At this time, output 92 is set to transfer control of 
the X stepper from mirror driver 70 to program 
controlled output 94 via gates 120 and 118. The X 
stepper is then ramped down in velocity to a stop 
in a number of pulses or steps equal to the fixed 

40 ramp-up number, e.g., 100 steps. The digital data 
in the two buffers may be processed at the end of 
a complete scan strip, as will be described below, 
or it may be processed as it is being digitized. At 
the conclusion of this scan strip, stepper motor 47 

45 moves the stage in the y direction so that a new 
scan strip can be run. Output 100 is used to send 
pulses to input 117a to step the motor; output 102 
determines whether movement is in the positive or 
negative y direction and output 104 determines the 

50 size of the steps (5 or 10u) and passes a signal to 
input 117c. Again, smaller or larger movements 
may be used under Protocol control. 

After stepping the Y motor 47 to move the 
specimen "up" or "down" in the y direction a 

55 distance equal to 60% of the scan length L, as 
shown in Fig. 4, the procedure described above, in 
which the X stepper motor 46 is ramped up in rate 
and moved the X distance, and ramped back down 
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to a stop is repeated but the stage is moved back 
to the left in the x direction. The Y stepper then 
moves the stage in the y direction a number of 
steps and the complete cycle is repeated. The 
number M of Y steps are counted by the program 
and the test is terminated when the Y distance (a 
given M) is reached as determined by the user 
through the appropriate Protocol parameter. 



System Setup and Calibration 



Upon starting the program, the user is asked to 
select a Protocol either from a set of named Pro- 
tocols stored on disk memory, which can be modi- 
fied at that time, or to generate a new Protocol. 
The Protocol is a set of parameters and demogra- 
phic information about the tests to be performed on 
the sample. Demographic information describes the 
type of sample, the number or name of the sample, 
and any comments relating to the sample. The 
Protocol, along with the cell count for the test, is 
stored with each set of test results to properly 
identify the data and test conditions. The Protocol 
parameters are: 1) the four start and end coordi- 
nates of the scan area and the X and Y step size, 
2) the analog data acquisition rate, 3) the number 
of sample points processed along each scan line, 
4) for each parameter, its name, its gain, a thresh- 
old level for cell findings a morphology threshold 
used to generate each cell's morphology word 
(binary neighborhood pattern), a flag to indicate 
which parameter will be used to find the cell's 
center, and where its values are displayed on the 
monitor during the test, 5) the size of the neighbor- 
hood (neighborhood is a concept to be described 
below), 6) the maximum cell count at which the test 
may be terminated, 7) the number of complete 
scans and 8) the prefix name of the file in which 
the test data is stored. 

Before beginning any sequence of tests, the 
user may be required to develop a set of calibra- 
tion data by methods described in further detail 
below. The user can select calibration data stored 
from a previous calibration run or can generate new 
calibration data. If new calibration data is desired, 
the user selects a section of the test specimen 
without cells or uses a specimen with no cells, 
places it on the microscope stage, and observes 
that the area to be scanned contains no cells or 
debris. The output of each parameter over the scan 
length is displayed continuously as a function of 
time on the monitor screen during calibration and 
stage position can be manually controlled. To de- 
termine these background values, the user keys 
the beginning of a background calibration scan and 
the values of 10 the parameters are digitized at one 
Y stepper position and one or more X stepper 



positions for a fixed number of scans on the order 
of 100. In one embodiment, this is accomplished in 
that the stage moves only a few steps in the x 
direction while continuously scanning up and down 

5 over generally the same blank (no sample) area. 
The parameter values of each of the scan sample 
point positions are averaged over all scans. The 
average values and their coefficients of variation 
are displayed on the monitor screen to be checked 

w by the operator to determine whether the machine 
is operating within the proper range of values. The 
coefficient of variation is the standard deviation of a 
measurement divided by the mean of that mea- 
surement. These values are stored in memory as 

75 an array to be used in calculating cell parameter 
values if required. 

The user can also initiate a second calibration 
run in which a test specimen provides an optical 
signal such as fluorescence or scatter. For exam- 

20 pie, a stained ground surface slide may be used. 
This is independently scanned as described above 
to derive a set of normalization factors for one or 
more parameters for each of the sample points. 
The calibration slide is scanned twice, the second 

25 time with the beam source blocked, so that the 
normalization factors are equal to the difference in 
the optical signal with and without the illuminating 
beam. This eliminates D.C. and other signal noise 
in the system that is present even in the absence 

30 of illumination. If a given parameter is not normal- 
ized its factor is one. These normalization factors 
are used to correct the data values for variations in 
the optical signal received from the cells due to, 
e.g., the cell's location on the slide, the variation in 

35 the incident laser light, or other variations not re- 
lated to a cell's constituents. These factors are then 
multiplied by the data points to remove the effect 
of these variations from the final measurements. 
The normalization factor may also be multiplied 

40 by a velocity normalization factor to account for the 
different velocities along each scan because if 
most of the scan is utilized, more samples are 
digitized for cells near the ends of each scan 
length L because the beam slows, stops and turns 

45 in the opposite direction before starting the next 
scan length L The velocity calibration factor is the 
cosine of the distance of a sample point from the 
center of a scan length and is multiplied by the 
normalization factor for each sample point and 

50 stored as an array to be used later to normalize 
each sample point value in calculating cell param- 
eter scan values. 

In another embodiment, an initial calibration of 
the device may be avoided by a method of self- 

55 calibration that operates during the sample run. 
This method avoids the need to scan a portion of a 
test specimen without cells or the use of any 
special calibration slide. The method generates an 
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illumination normalization factor used to calibrate 
the optical data to correct for beam illumination and 
light collection differences along each scan line 
(length L) across the sample. When the cells of the 
sample are used to generate the normalization fac- 
tor, the method is self- calibrating. In the alter- 
native, however, the method can be used with 
calibration particles in an initial run in which the 
calibration particles are scanned for a constant 
optical value of one of a number of known constant 
parameters, such as size or flurescence, to gen- 
erate the illumination normalization factor. 

This method may be used when the cell popu- 
lation of the sample includes a number of cells that 
have a constant optical value for some parameter, 
such as, for example, the io DNA value of cells in 
the resting state or the size of the cells. To achieve 
an accuracy of about three percent in the cali- 
brated data, approximately one thousand cells of 
the population having this constant parameter value 
should be scanned and counted. 

For example, in a typical population of cells, a 
large fraction of those cells are in the resting state 
(G ) in which all of the cells have the same DNA 
value. In the ideal population measurement, a sub- 
stantial number of cells would all have the same 
value for this particular parameter. Using this ideal 
as a goal, the optical data is calibrated by the 
illumination normalization factor to obtain a new set 
of data with a distribution as close to this ideal 
result as possible. A normalization factor can be 
developed and applied to the present sample, to 
subsequent samples, or, each new sample run can 
be corrected based on a normalization factor devel- 
oped from all prior runs. 

In particular, this calibration method functions 
as follows (as illustrated in Fig. 9): As the cells of 
the sample are scanned, a data list is generated for 
each cell including, inter alia , the optical neighbor- 
hood value (step 300) and the location of the 
maximum value for a select parameter (step 302), 
which represents the center of the cell along the 
scan line, for one or more select parameters. The 
neighborhood is then recentered with respect to 
this maximum value (step 304). In this embodi- 
ment, the optical neighborhood value is not cali- 
brated for other factors before being multiplied by 
the illuminization normalization factor. Next, the op- 
tical neighborhood value occurring with the highest 
frequency for the complete sample run is deter- 
mined and designated as the "mode optical value" 
(step 306). Next, the optical neighborhood values 
are reviewed and ail cells having optical neighbor- 
hood values within a predetermined range of plus 
or minus "D" of the mode optical value are se- 
lected to generate a selected subpopuiation of cells 
(step 308). This subpopuiation of cells is then sort- 
ed as a function of location, or pixel position, of the 



maximum value, or center, of each cell along the 
scan line (step 310). The value of D is user defined 
and is typically set to between 2 and 6 percent of 
the mode optical value. If this range is selected to 

5 be too broad, it will select cells that do not have 
the constant value, and if this range is too small, 
then not enough cells to create an accurate aver- 
age will be selected. When calibration particles that 
are known to be uniform are used, there is no need 

w to select a subpopuiation of particles within a pre- 
determined range around the mode optical calibra- 
tion value, because ail the particles are the same. 

Next, the average optical value of the selected 
subpopuiation of cells for each pixel position along 

75 the scan line is determined (step 312). Thereafter, 
an array of illumination correction coefficients is 
computed for each scan line including one correc- 
tion coefficient for each location, or pixel position 
along a scan line (e.g., for 18 pixels that are 

20 preferably used along a scan length L of 30 sam- 
ples, there would be 18 coefficients) (step 314). 
Each correction coefficient in the array equals the 
ratio of the corresponding mode optical value to 
the average optical value for the associated param- 

25 eter. This array is the illumination normalization 
factor. 

The array of correction coefficients is one di- 
mensional in that it corrects for variations in light 
excitation or collection along the scan line. How- 

30 ever, the method may also be used to generate 
such an array for two-dimensional systems. For 
example, the laser beam may be incrementally 
moved after each scan over some range in a 
direction perpendicular to the present scan beam, 

35 rather than moving the microscope stage. At the 
end of this range, the beam movement would be 
reset to zero and the stage incremented. In this 
case, the resulting optical data would be corrected 
as determined by peak location, i.e., the center of 

40 the cell in two dimensions. 

All summed values, i.e., the optical neighbor- 
hood values, for each parameter are then calibrated 
for differences in illumination along each scan line 
by finding the peak pixel position for each cell and 

45 multiplying the optical neighborhood value by the 
normalization factor for that scan line position for 
each parameter (step 316). This generates a cali- 
brated set of optical data. 

It is desirable to repeat this process for mul- 

50 tiple iterations until the data stabilizes. This is ac- 
complished by using the calibrated optical data 
generated by the first run as the basis for second 
and subsequent calibrations. In other words, the 
calibration method is applied to the initially cali- 

55 brated optical data stored in the computer memory 
to generate a new and more accurate illumination 
normalization factor, which is then multiplied by the 
initially calibrated optical data to generate a second 
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set of more accurately calibrated optical data. This 
process can be repeated many times, although the 
data should stabilize after approximately three iter- 
ations. The illumination normalization factor can be 
applied directly to subsequent sample run optical 
data or the factor can be updated by applying the 
method to each run and averaging the latest set of 
correction coefficients from the coefficients of prior 
runs and applying the averaged set to the next run. 



Optical Signal Data Accuisition 

After selecting a Protocol and calibrating the 
instrument, the operator may initiate one or more 
data acquisition runs. The purpose of these runs is 
to scan an area of a test slide defined in the 
Protocol, find all of the cells on the slide that meet 
a threshold criteria for a given parameter that is set 
in the Protocol by the operator, and generate a list 
for each cell found that contains the corrected 
optical and morphology values for each parameter 
of the Protocol, the measurement time, and the X 
and Y position of each cell with respect to Home. 
The user places the test specimen on the stage 
and initiates a run by typing a key on the computer 
keyboard. The flow chart of Fig. 3 illustrates the 
general mechanical optical signal data acquisition 
loop as described above. 

The program causes the stage to be driven to 
Home and then moves the stage over the test area. 
As the stage moves, the beam is scanned back 
and forth to create a scan path as shown in Fig. 4. 
For each X direction scan strip, the optical signal 
value of each parameter provided by the specimen 
is digitized by an analog-to-digital (A/D) converter 
at a sample rate set by the Protocol to create a 
sequence of digital data sample values. Typically 
120,000 such optical signal values are digitized and 
stored in a buffer memory within 3 seconds. Typi- 
cally each scan contains 20 samples so that 2000 
to 6000 forward and reverse scans corresponding 
to two to four parameters with steps of 5 to 10 
microns cover a length of the same order of mag- 
nitude as a slide coverslip. Data acquisition and 
processing may be done sequentially or simulta- 
neously. 

The Relationship Between Spot Size and Sampling 
RaTe 

In the described embodiment, the optical sys- 
tem is designed so that the spot size of the scan- 
ning laser beam spot 12c at the specimen surface 
28 is approximately 10 microns when a 20X micro- 
scope objective is used in the revolving nosepiece. 
The spot size can be changed by rotating the 



nosepiece to 40X to produce a 5 micron spot or to 
10X to produce a 20 micron spot. 

As the beam 12b scans across surface 28, a 
continuous optical analog data stream is produced 

5 by each of the instrument's optical sensors 20, 22, 
24, 26. Figs. 5a through 5d show graphically how 
the analog optical data frm a fluorescence sensor 
is sampled by an A/D converter. Fig. 5a shows a 
hypothetical profile of a fluorescence signal from a 

m cell with a level C representing a base level of a 
constituent and with levels representing constituent 
variations Ci and C2. Fig. 5b shows a scan spot 
intensity profile. The "size" of the spot is measured 
at the half-power points of the Gaussian profile. 

75 Fig. 5c shows a digital representation of the flu- 
orescence signal with inadequate sampling, and 
Fig. 5d shows a digital representation of the flu- 
orescence signal with sampling according to the 
invention. The arrows represent the sampling inter- 

20 val. 

The user can define the rate at which each of 
these continuous optical analog data streams are 
sampled and converted by the A/D converter and 
stored as digital data values in the computer mem- 

25 ory. The rate can be set using the hardware pro- 
tocol up to a maximum rate of 100,000 samples 
per second. A typical experiment uses a rate of 
32,000 dual parameter samples per second, i.e., 
the two analog data streams are each sampled at 

30 31 .25 microseconds. 

The velocity of scan beam 1 2b is fixed by the 
resonant frequency of the resonant scanner 70. 
Such scanners are commercially available with fre- 
quencies to 2400 Hz. One embodiment uses an 

35 800 Hz. scanner, i.e., scan beam 12b moves for- 
ward and back, each with a distance "L", the scan 
length, 800 times per second, taking 625 micro- 
seconds for the beam to move from one end to the 
other of the scan length L with a cosine velocity 

40 profile. (625usec/31 .25usec = 20 the number of 
sample points per scan length.) The scan path of 
the laser beam is shown in Fig. 4. 

As shown in Fig. 5b the intensity profile of the 
scan spot 12b on the specimen surface 28 is 

45 Gaussian. The analog data generated by any sen- 
sor from a continuous scan of an optical signal 
from the specimen is the convolution of the Gaus- 
sian intensity profile and the distribution of the 
desired cellular constituent signal, i.e., the two dis- 

50 tributions are multiplied as one is shifted in time. 

Consider a case in which a fluorescence dyed 
cell is scanned. A hypothetical profile of fluores- 
cence is shown in Fig. 5a which represents levels 
of fluorescence C due to the cell's base level of a 

55 constituent and fluorescence due to constituent 
variations C1 and C 2 . If the sampling rate is low so 
that the beam moves more than a spot width 
between the locations at which the analog data 
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stream is sampled, some of the information in the 
fluorescence signal will be lost. The digital samples 
will be added later to try to determine the amount 
of dye along a line segment. In the case of low 
sampling, i.e., samples are taken too infrequently, 5 
the resultant set of values would not be representa- 
tive of the convolution of the distributions of the 
beam intensity profile and the dye distribution be- 
cause the sampling rate is too low producing 
"gaps" in the digital data corresponding to the 10 
analog data as shown in Fig. 5c. Summing these 
samples would underrepresent parts of the base 
level C and furthermore, the information represent- 
ing the actual increase due to constituent variation 
Cz is lost between the two peaks representing base 15 
level fluorescence C. If the sampling rate is such 
that successive samples in the specimen plane are 
spaced apart by the spot size or less, all parts of 
the specimen will be equally represented as shown 
in Fig. 5d, because Gaussian profile spots spaced 20 
one spot size apart when added produce a uniform 
distribution along the scan. This is because the 
Gaussian profiles overlap enough when spaced one 
"spot size" apart that, when added, the total equals 
the maximum or peak value of the Gaussian curve. 25 
The cell specimen is sampled with close enough 
spot spacing so that when the samples from a 
fluorescence sensor are added, with proper calibra- 
tion, it will yield the accurate amount of. the cells' 
constituents independent of the distribution of that 30 
constituent in the cell or the location of the cell on 
the slide. 

The described embodiment uses a conven- 
tional circuit card available from Data Translation 
(DT) for converting analog data to digital values. 35 
This card utilizes direct memory access (dma) in 
which the DT card is set up by the program to 
continuously sample each of up to four analog 
signals at a predetermined rate. This DT card also 
converts each of these analog samples to its digital 40 
data value equivalent, and to store each set of 
signal samples in successive memory locations 
with a predetermined starting memory location, 
completing the operation when a predetermined 
number of memory locations (buffer) is filled with 45 
digital values. 

The program sets up this analog-to-digital op- 
eration by specifying the signals to be sampled, 
the sampling rate, the initial memory location, and 
the size of the buffer. One complete X strip is 50 
processed from raw data into list data in one batch 
so that the buffer size must be equal to the number 
of scans in a strip (typically 2500) times the num- 
ber of samples in a scan (typically 20) times the 
number of optical signals used in the experiment 55 
(typically 2 to 4); that is, the buffer size is approxi- 
mately 100,000. When the buffer is full, program 
control is transferred to the functions that will find 



cells, find their peak values, and determine their 
optical and morphological property values. At this 
point, there is a buffer with typically on the order of 
100,000 sequential values in it. The dma operation 
and list value determination operations can also be 
overlapped. 



Digital Data Processing 

Once the digital data sample values are stored 
in the computer memory (buffer) the data is pro- 
cessed by a variety of protocol controlled functions 
to correct the data, e.g., for background, to cali- 
brate the data, and to generate the desired optical 
and morphological property values. These function 
steps are shown in the flow chart of Fig. 7. 

The first data processing program function 
step, 200 in Fig. 7, locates the beginning of each 
forward scan (typically 2500 per strip) in.the strip. 
As shown in Fig. 6, one of the signals, typically the 
forward scatter signal, has added to it a synchro- 
nization pulse of sufficient amplitude, duration, and 
negative polarity so that it can be distinguished 
from all normal signals. The synchronization pulse 
is derived directly from the mirror scan driver 70, 
occurring once for each forward scan at the same 
time for every scan and near its beginning so that 
this pulse does not interfere with the actual data. In 
the described embodiment, only the middle 
(approximately 12) sample points of the 20 sample 
scan (i.e., there are 20 samples per scan length L) 
are used and will be referred to as POINTS. A 30 
sample scan length Would have approximately 18 
POINTS. The program first searches successive 
memory locations of the buffer for sequences of 
three values that match the shape of the synchro- 
nization pulse and produces a table of pointers to 
all of these locations. The starting location of every 
data POINT value is fixed at a known displacement 
in the buffer from each of these pointer locations. 
Because the relationship of the synchronization sig- 
nal and scan position is fixed, the start position of 
every scan can be marked and the data buffer 
values can each be associated with a specific scan 
position by appropriate record keeping in the pro- 
gram. 

Groups of contiguous points for each param- 
eter are located in memory based upon their dis- 
tance from the synchronization signals. If param- 
eter samples are taken sequentially, i.e., POINT 1 
for parameter A, POINT 1 for parameter B, POINT 
2 for parameter A, etc., then these groups of con- 
tiguous points will include every other sample in 
memory (for 2 parameters). 

The next function step, 202 in Fig. 7, called by 
the program is responsible for finding cells. The 
event of locating an individual cell is called a "cell 
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find". The function causes a search of each of the 
POINTS of each scan to find values for any param- 
eter (with its background value subtracted) that 
exceed the Protocol threshold value for that param- 
eter. Thus, any one parameter can be used to find 
a cell by an increase in that given parameters 
signal value caused by that cell. 

Background values can be obtained from a 
calibration area as described above or they can be 
generated, either by averaging or determining the 
minimum values, of samples in the same position 
of scans in one or more prior or subsequent neigh- 
borhoods, initially, the program computes a table of 
initial background values for the POINTS of a part 
of a scan strip, i.e., the first approximately 50 
scans, to determine the lowest value for ail 
POINTS for each parameter signal. The initial back- 
ground value for these POINTS is stored in a table. 

One of the Protocol parameters of this inven- 
tion is the size of the neighborhood. The neighbor- 
hood has a value N that is typically near five 
POINTS but is adjusted by the user to the size and 
concentration of cells in the test. The neighborhood 
is set to statistically include only one cell. 

In the cell find function, a separate step 204 in 
Fig. 7, defines a neighborhood around each cell 
found. The neighborhood value N defines the area 
(POINTS x scans) surrounding the cell find in 
which other cell finds are prohibited, that is, once a 
single cell is found in a given neighborhood, the 
next find must be outside that neighborhood. The 
number of scans skipped after the end of a neigh- 
borhood is typically 3. Also, because data values 
are added only within the neighborhoods, cell finds 
are not allowed in areas of the sample at the scan 
edges which comprise less than a full neighbor- 
hood. The present real time, including date, is 
taken from the computer clock and stored in the 
cell list. The value of N, which is equal to a number 
of sample POINTS of a neighborhood and the 
number of successive scans of a neighborhood is 
defined in the Protocol. N is preferably odd for 
ease of defining a center of the neighborhood. 

The next function step, 206, compares all val- 
ues along N POINTS of N forward scans, i.e., 
within one neighborhood, centered at the location 
of the cell find to determine the maximum value of 
that POINT minus its corresponding POINT back- 
ground value for one of the Protocol determined 
signal parameters. The maximum value defines the 
center of the cell and the new location of the cell 
find. The neighborhood is recentered in step 207. 
The number of the scan, the POINT number, and 
the strip number are all stored in the main cell list 
file to define this new cell location, along with the 
computer system time. 

In the next step, 208, a new set of POINT 
background values for each neighborhood is com- 



puted by determining the minimum values of sam- 
ples for each parameter signal using a set number 
of scans (typically 5) that end outside the neighbor- 
hood, e.g., just prior to the beginning of the current 

5 neighborhood. Subtracting this background from 
the digital data POINTS generates optical data 
POINTS in the computer memory. 

Intensity, velocity and background as deter- 
mined by the calibration run and by the program 

70 are all taken into account in the next step 210 by 
multiplying each POINT value for each parameter 
by an array of POINT intensity correction factors 
and velocity correction factors. The next function 
step, 212 in Fig. 7, sums all corrected optical data 

15 values for each parameter for each sample of the 
defined neighborhood centered at the position of 
the peak value for that cell. This sum is the optical 
neighborhood value. In the self-calibration method 
described above, this optical neighborhood value is 

20 not initially calibrated for intensity or velocity, but is 
corrected in a subsequent step. 

In this embodiment, the number of samples of 
the cell taken in each direction is greater than the 
projection of cell length in that dimension divided 

25 by the spot size in that direction so that when 
these samples are added, the result accurately 
represents the integrated value of that cell param- 
eter and can be directly related to the absolute 
quantity of a given cell constituent. This results 

30 from the principle that the samples are taken at a 
distance apart that is equal to or less than the spot 
size of the incident light beam. The spot size is 
selected by choice of microscope objective power, 
whereas the sampling rate, the step size, and the 

35 size of the neighborhood are all selected in the 
Protocol to adjust the test to the speed, accuracy 
and sizes of cells tested in the specific test. 

The program in the next step, 214, generates a 
digitized binary neighborhood pattern or morphol- 

40 ogy word for each cell describing its shape. Each 
POINT in the neighborhood for each parameter is 
compared with that parameter's "morphology 
threshold" set in the Protocol. For each POINT 
where the value exceeds the morphology thresh- 

45 old, a specific position of the morphology word in 
memory is set from zero to one and the word is 
shifted by one. If the threshold is not exceeded, the 
location is kept at zero and the word is shifted by 
one. Thus, each parameter of each cell is used to 

so generate a word equal in length to the area of the 
neighborhood (POINTS x scans) that contains a 
digital image of the cell. In practice, the size of this 
word can be made smaller than the number of 
POINTS to reduce storage space. In a final step, s 

55 216, these words are stored in a list for each cell 
along With data for the neighborhood value gen- 
erated in step 210, the time of the scan, and cell 
position data. 



13 



25 



EP 0 421 736 A2 



26 



Display of Optical and Morphology Parameters 

The above set of functions is repeated for each 
cell find in the buffer representing one X strip. The 
optical property values (generated from the optical 
neighborhood value) or the morphological property 
values (generated from the binary neighborhood 
pattern) of two parameters are displayed on a 
monitor screen as a dot whose x position is propor- 
tional to one parameter and whose y position is 
proportional to a second parameter. When a spe- 
cific number K of cells are found (typically 500) or 
a specific number of Y-steps (M) is reached, the 
"run" is completed and a list of values for each cell 
is written in a file to a fixed disk of the computer. 
The device also may be programmed to repeat a 
run, or rescan, for any number P of passes. 

To repeat a run, the microscope stage is 
moved to begin scanning the first X strip and all of 
the above steps are repeated. Simultaneously, oth- 
er program functions control the position of the 
stage. As in the first pass, the stage is initially 
moved to a HOME position, then to an initial start 
position and ramped up to scanning speed while 
keeping track of position, moved in synchronized 
steps across a strip in the x direction, ramped 
down in speed to a stop, and stepped in the y 
direction as described above. In this way, by re- 
peat scanning of the same specimen of cells, the 
operator can observe kinetic changes in the individ- 
ual cells. 

The list of data generated as described above 
can be processed by a data display program that 
can be run at a later time or between strips of the 
above data acquisition program. The principal fea- 
tures of this program are as follows. 

Using a "Data Protocol" the user can define a 
list of properties that are computed by applying 
mathematical functions to one or two optical neigh- 
borhood values, morphology words or time param- 
eters. For example, an optical property can be an 
optical neighborhood value for one parameter (e.g., 
a value that can be converted into a number of 
picograms of DNA), the ratio of two such values, 
the time, or the sum of all of the 1 bits of a 
morphology word, i.e., binary neighborhood pattern 
(the cell area). These properties are used as the 
basis for a sequence of cell selections to count and 
display cell subpopulations on the monitor screen. 

Screens are generated on the monitor with 
axes corresponding to each of two selected prop- 
erties. Each cell is displayed on the screen as a 
dot whose position on the axis is proportional to a 
property and whose color is proportional to the 
number of cells occurring at that screen position. 
Using a user controlled mouse, for example, a user 
can bound a selected area (this can be extended to 
additional areas) of the screen by a polygon. Two 



additional properties that may be the same or 
different from those of the first screen are des- 
ignated in the data protocol. Cells whose properties 
are within the area of the polygon of the first 

5 screen are counted and displayed on a second 
screen whose coordinates are the second set of 
properties. A number of polygons, e.g., six, can be 
drawn on the second screen and cells within each 
of these regions are counted. As shown in Fig. 8, 

w the counts appear on the screen next to their 
corresponding regions. 

In addition to counting subpopulations of cells 
based on successive counts of derived properties, 
the user can display selected data in a variety of 

75 other formats including isometric projections or 
contour maps of any set of properties. Another 
aspect of the instrument according to the invention, 
as illustrated in Fig. 8, is the ability of this program 
to generate a polygon on the monitor screen that 

20 represents a specific location or set of parameter 
values. The program can then instruct the instru- 
ment to locate each cell whose properties fall with- 
in the polygon area. That is, each polygon contains 
a subpopulation of cells within the sample that 

25 exhibits a specific characteristic common to each 
of the cells within the polygon. As each cell is 
automatically positioned under a microscope objec- 
tive 39 for visual observation, a dot on the screen 
corresponding to that ceil within the population is 

30 marked by a blinking circle centered on it (as 
shown in the lower right hand polygon in Fig. 8). 

After completion of each X strip, the values for 
each cell's parameters are added to a monitor 
screen display. During each test, the user sees 

35 each cell represented as a dot on the screen at a 
position determined by the values of two param- 
eters or as two such patterns if more than two 
parameters are used. The axes of these graphs are 
labeled using the Protocol parameter names. Ac- 

40 cording to the invention, the positions of found cells 
may also be represented by dots on an additional 
graph on the screen. After a fixed number of cells 
is found, typically about 500, the data lists for 
those cells are stored in a disk file containing a 

45 header with the Protocol. After either the total area 
has been scanned or the Protocol maximum cell 
count is reached, the remaining cell data is stored 
in a list file followed by the total cell count. 

The user can call for the program to repeat the 

so scan one or more times. In these subsequent 
passes, all of the above functions, with the optional 
exception of cell finding, are repeated. The loca- 
tions of all found cells from the first pass are read 
into a table in memory and these locations are 

55 used at the appropriate times to indicate the start- 
ing position of the peak value finding routine. The 
new time is stored in the new list and all other 
functions are identical to the first pass. 
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Data Analysis 

The result of the program of the invention is a 
file containing a set of data for typically 1000 to 
100,000 cells. These data values are used as the 
input to an analysis and display program which can 
be identical to any of the programs described in 
the flow cytometry literature. These programs are 
reviewed in the texts, Melamed, et al. Flow 
Cytometry and Sorting (John Wiley & Sons, 1979, 
Ch. 20), or~SFapiro, H., Practical Flow Cytometry . 
(A.K. Liss, 1985, Ch. 5). 

The monitor screen displays of the present 
invention describing morphology, time and position, 
however, are another feature of the present inven- 
tion. For example, time can be represented graphi- 
cally as one additional parameter along with scatter 
or fluorescence parameters or specific functions of 
parameters may be tested against time variable 
thresholds or displayed as a function of time to 
implement experiments or tests. Time can also be 
displayed by modifying dot or contour colors on a 
display to mark the passage of time. Functions of 
the morphology word such as the total one bit 
count indicative of area of that parameter may be 
used like any parameter. It is contemplated that the 
large number of total parameters that can be ob- 
tained from the morphology and value data can be 
best analyzed with techniques described in the 
pattern recognition or artificial intelligence literature 
to achieve the end goal of identifying specific ceil 
types or pathological cells on a slide. 

Once cells with unique parameter values are 
found by the analysis program, the test specimen 
can be placed on the microscope stage and a 
program called to automatically center the stage on 
each of these unique cells in turn by reading in a 
list of the X and Y values of cells to be reviewed. 
Properties of these cells may be displayed on the 
monitor screen. The user may use this to verify the 
automatic data analysis, examine the cells to fur- 
ther determine their pathology, check the quality of 
the manual procedures in a clinical laboratory, or 
physically separate or sort cells of interest using 
cookie cutter like devices. 

The kinetic behavior of live cells, primarily 
blood lymphocytes, when reacted with certain anti- 
gens during a test, may be an important new 
method of testing for systemic pathologies such as 
autoimmune disease or cancer. As described 
above, cells with specific properties may be iden- 
tified and their locations stored in memory. The 
specimen in a cuvette or chamber on the stage can 
be rescanned and the new values of these selected 
cells recalculated and listed. 

The device according to the invention also con- 
templates sequential analysis techniques to maxi- 
mize speed of analysis. Cells of interest found 



during a first pass will be examined for other pa- 
rameters or at greater resolutions to determine 
properties that require greater analysis complexity 
and time. 

5 Other embodiments are feasible. 



Claims 

10 1. A method for generating optical data that char- 
acterizes a population of cells, the method com- 
prising: 

(a) scanning the ceil population with a beam to 
produce digital data samples, different said sam- 

75 pies representing different locations within the 
cell population; 

(b) storing said digital data; 

(c) locating a cell within said population; 

(d) defining a neighborhood around the located 
20 cell; 

(e) estimating a background level for the neigh- 
borhood based upon stored digital data cor- 
responding to locations outside said neighbor- 
hood; and 

25 (f) correcting each of said digital data samples 
corresponding to the neighborhood with the es- 
timated neighborhood background level to gen- 
erate the optical data, said sample 

2. The method of claim 1, wherein said sample 
30 correcting step comprises subtracting the corre- 
sponding background level from each of said 
stored digital data samples in a neighborhood. 

3. The method of claim 1 , wherein said digital data 
samples represent a plurality of morphological 

35 properties of each cell in the population. 

4. The method of claim 1, wherein said neighbor- 
hood comprises a given number of sample points. 

5. The method of claim 4, wherein said background 
level for a neighborhood comprises a plurality of 

40 individual background levels, one for each sample 
point in said neighborhood. 

6. The method of claim 1 wherein said ceil popula- 
tion is on a movable surface, said scanning pro- 
ceeds in synchrony with a stepping motion of said 

45 surface and digital location data is generated from 
which the location corresponding to a sample can 
be derived. 

7. The method of claim 6, wherein said cells on 
said movable surface can be visually observed. 

50 8. The method of claim 1, wherein said neighbor- 
hood background level is determined by finding 
minimum digital data sample values in memory 
locations corresponding to actual locations adjacent 
the neighborhood for each sample in said neigh- 

55 borhood. 

9. The method of claim 1, wherein said neighbor- 
hood background levels are determined by averag- 
ing a plurality of digital data sample values in 
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memory outside said neighborhood. 

10. The method of claim 1, wherein said cell popu- 
lation is on a movable surface, said scanning pro- 
ceeds in synchrony with a stepping motion of said 
surface, the time corresponding approximately to 
the time the digital data was acquired is recorded 
and stored as synchronous digital time data, and 
digital location data is generated from which the 
location corresponding to a sample can be derived, 
whereby each sample has a corresponding digital 
time point and location. 

11. The method of olaim 10, further oomprising the 
step of rescanning the population of cells, generat- 
ing new digital data, location data and time data, 
and comparing said new data with the correspond- 
ing data from a previous scan to determine 
changes in said cell population. 

12. The method of claim 1, wherein the size of said 
neighborhood is set to statistically include a single 
cell. 

13. The method of claim 1, wherein said locating 
step comprises detecting a digital data sample 
corresponding to the edge of the cell first con- 
tacted by said beam. 

14. A method for generating optical data that char- 
acterizes a population of cells, the method com- 
prising: 

(a) scanning the cell population with a beam 
with a spot size comparable in size to the cells 
to be scanned to produce digital data samples, 
different said samples representing different lo- 
cations within the cell population and said sam- 
ples being taken at a sampling rate set such that 
the distance travelled per time as the spot is 
scanned between consecutive sampled locations 
approximates said spot size; and 

(b) processing the digital data to generate the 
optical data that characterizes the cells in the 
population 

15. The method of claim 14, wherein said spot size 
is 3 to 50 microns. 

16. The method of claim 14, wherein said process- 
ing step comprises correcting said digital data 
samples with a background level to generate said 
optical data. 

17. A method for generating calibrated optical data 
that characterizes a population of cells, the method 
comprising: 

(a) scanning the cell population on a surface 
with a beam to produce digital data samples, 
said beam having a given illumination intensity, 
different said samples representing different lo- 
cations within the cell population; 

(b) storing said digital data; 

(c) scanning a region of said surface without 
cells to generate intensity calibration data; 

(d) generating an intensity normalization factor 
from said intensity calibration data and storing 



said intensity normalization factor in memory; 

(e) locating a cell within said population; 

(f) defining a neighborhood around the located 
cell; said neighborhood containing optical data 

5 derived from said stored digital data samples 
within said neighborhood; and 

(g) calibrating the optical data in said neighbor- 
hood with said intensity normalization factor to 
generate the calibrated optical data. 

w 18. The method of claim 17, wherein said optical 
data calibrating step comprises multiplying said 
optical data by said intensity normalization factor. 

19. The method of claim 17, wherein said surface 
without cells is a calibration slide comprising a 

15 uniform fluorescent dye surface. 

20. The method of claim 17, wherein said surface 
without ceils is a calibration slide comprising a 
uniform scattering surface. 

21. The method of claim 17, wherein the derivation 
20 of said optical data contained in said neighborhood 

comprises the steps of: 

estimating background levels for the individual 
neighborhoods based upon stored digital data cor- 
responding to locations outside said nefghborhood; 
25 and correcting each of said digital data samples in 
each neighborhood with the corresponding esti- 
mated neighborhood background level to generate 
said optical data. 

22. A method of measuring cellular properties com- 
30 prising the steps of: 

(a) scanning the cell population with a beam to 
produce digital data samples, different said sam- 
ples representing different locations within the 
cell population; 
35 (b) storing said digital data; 

(c) locating a cell within said population; 

(d) defining a neighborhood around the located 
cell, said neighborhood containing optical data 
derived from said stored digital data samples 

40 within said neighborhood; 

(e) summing all samples of optical data in a 14 
neighborhood to generate an optical neighbor- 
hood value for that neighborhood; 

(f) comparing each sample in a neighborhood 
45 with a preselected morphology threshold to gen- 
erate a binary neighborhood pattern of samples 
either above or below said threshold; 

(g) storing said optical neighborhood value and 
binary neighborhood pattern in memory; and 

50 (h) generating optical properties corresponding 
to said optical neighborhood value and mor- 
phological properties corresponding to said bi- 
nary neighborhood pattern for each cell. 

23. The method of claim 22, wherein the genera- 
55 tion of said optical properties comprises: applying 

math functions to one or more said optical neigh- 
borhood values to generate optical properties of 
the cell corresponding to said neighborhood. 
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24. The method of claim 22, wherein the genera- 
tion of said morphological properties comprises: 
applying boolean functions to one or more said 
binary neighborhood patterns to generate morpho- 
logical properties of the cell corresponding to said 
neighborhood. 

25. The method of claim 22, wherein the derivation 
of said optical data contained in said neighborhood 
comprises: 

estimating a background level for the neighborhood 
based upon stored digital data corresponding to 
locations outside said neighborhood; and 
correcting each of said digital data samples cor- 
responding to the neighborhood with the estimated 
neighborhood background level to generate said 
optical data. 

26. The method of any one of claims 1, 14, 17, or 
22, wherein said digital data samples produced by 
the cells 3 and the beam are derived from analog 
optical data. 

27. The method of claim 26, wherein said optical 
data is fluorescence, forward angle light scatter, 
extinction, or wide angle light scatter. 

28. An apparatus for measuring cellular properties 
of a population of cells comprising: 

a beam source for generating a beam of elec- 
tromagnetic radiation to illuminate the cells to be 
measured with a beam spot comparable in size to 
the cells to be scanned, an optical signal being 
generated as a result of said spot illuminating said 
cells; 

a surface upon which the cells are located; 
an optical path for directing said spot from said 
beam source to the cells on said surface; 
a scanner interposed in said optical path between 
said beam source and the cells for scanning said 
spot across said surface so that said spot passes 
over a plurality of cells on said surface; 
a sensor for measuring said optical signal; 
an analog-to-digital converter arranged to sample 
said optical signal at a specified rate and to pro- 
duce digital data; whereby said sampling rate is set 
such that the distance traveled per time as the spot 
is scanned between consecutive sampled locations 
approximates said spot size; 
a data processor for deriving optical data from said 
digital data, generating optical property values from 
said optical data corresponding to individual cells, 
and processing said values to determine cellular 
properties; and a memory for storing said digital 
data produced by said converter. 

29. The apparatus of claim 28, wherein said spot 
size is 3 to 50 microns. 

30. The apparatus of claim 28, wherein said scan- 
ner is a resonant galvanometer scanner. 

31. The apparatus of claim 28, wherein the 
wavelength distribution of said beam is controlled 
by the data processor. 



32. An apparatus for generating optical data that 
characterizes a population of cells comprising: 

(a) a scanner for scanning the cell population 
with a beam to produce digital data samples, 

5 different said samples representing different lo- 
cations within the cell population; 

(b) a memory to store said digital data; 

(c) means for locating digital data samples in 
memory corresponding to a cell within said pop- 

70 ulation; 

(d) means for defining a neighborhood around 
the located cell; 

(e) means for estimating a background level for 
the neighborhood based upon stored digital data 

rs corresponding to locations outside said neigh- 
borhood; and 

(f) means for correcting each of said digital data 
samples corresponding to the neighborhood with 
the estimated neighborhood background level to 

20 generate the optical data. 

33. An apparatus for generating calibrated optical 
data that characterizes a population of cells com- 
prising: 

(a) a scanner for scanning the cell population on 
25 a surface with a beam to produce digital data 

samples, said beam having a given illumination 
intensity, different said samples representing dif- 
ferent locations within the cell population, and 
for scanning a region of said surface without 
30 cells to generate intensity calibration data; 

(b) a memory for storing said digital data; 

(c) means for generating an intensity normaliza- 
tion factor from said intensity calibration data; 

(d) means for locating digital data samples in 
35 memory corresponding to a cell within said pop- 
ulation; 

(e) means for defining a neighborhood around 
the located cell; said neighborhood containing 
optical data derived from said stored digital data 

40 samples within said neighborhood; and 

(g) means for correcting the optical data in said 
neighborhood with said intensity normalization 
factor to generate the calibrated optical data. 
34. An apparatus for measuring cellular properties 

45 comprising: 

(a) a scanner for scanning the cell population 
with a beam to produce digital data samples, 
different said samples representing different lo- 
cations within the cell population; 

50 (b) a memory for storing said digital data; 

(c) means for locating digital data samples in 
memory corresponding to a cell within said pop- 
ulation; 

(d) means for defining a neighborhood around 
55 the located cell; 

(e) an adder for summing all samples of optical 
data in a neighborhood to generate an optical 
neighborhood value for that neighborhood; 
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(f) a comparator for comparing each sample in a 
neighborhood with a preselected morphology 
threshold to generate a binary neighborhood 
pattern of samples either above or below said 
threshold; and 

(g) means for generating optical properties cor- 
responding to said optical neighborhood value 
and morphological properties corresponding to 
said binary neighborhood pattern for each cell. 

35. A method for generating calibrated optical data 
that characterizes a population of cells in which a 
number of the cells has a constant optical value, 
the method comprising: 

(a) scanning the cell population with a beam 
along a scan line to produce digital data sam- 
ples, different said samples representing dif- 
ferent locations along said scan line within the 
cell population; 

(b) storing said digital data; 

(c) locating a cell within said population; 

(d) defining a neighborhood around the located 
cell, said neighborhood containing optical data 
derived from said stored digital data samples 
within said neighborhood; 

(e) determining the maximum digital data sam- 
ple within said neighborhood and recentering 
the neighborhood with respect to said maximum 
digital data sample; 

(f) summing all optica! data values in a reen- 
tered neighborhood to generate an optical 
neighborhood value for that neighborhood; 

(g) processing said optical data representing 
individual cells to generate an illumination nor- 
malization factor to correct optical data for illu- 
mination differences as a function of cell loca- 
tion along the scan line; and 

(h) correcting optical data for each cell at dif- 
ferent locations along the scan line by multiply- 
ing the optical neighborhood value for that cell 
by the corresponding illumination normalization 
factor. 

36. The method of claim 35, wherein said process- 
ing of said optical data comprises: 

(a) determining the most frequently occurring 
optical neighborhood value to define a mode 
optical value; 

(b) selecting a subpopulation of cells with optical 
neighborhood values within a predetermined 
range around said mode optical value; 

(c) determining a set of average optical neigh- 
borhood values of cells in the subpopulation as 
a function of cell location along the scan line; 
and 

(d) computing an array of correction coefficients 
equal to the ratio of the mode optical value to 
the average optical value for each location along 
the scan line to generate an illumination nor- 
malization factor. 



37. The method of claim 36, wherein said predeter- 
mined range of optical neighborhood values select- 
ing said subpopulation of cells is from plus or 
minus 2 to 6 percent of said mode optical value. 
5 38. The method of claim 36, wherein said sub- 
population of cells comprises cells in the resting 
phase (Go). 

39. The method of claim 36, wherein said sub- 
population of cells comprises cells of substantially 

10 the same size. 

40. The method of claim 36, wherein the steps of 
computing an array of correction coefficients equal 
to the ratio of the mode optical value to the aver- 
age optical value for each location along the scan 

15 line to generate an illumination normalization factor, 
and correcting optical data for each cell at different 
locations along the scan line by multiplying the 
optical neighborhood value for that cell by the 
corresponding illumination normalization factor, are 

20 repeated iteratively for initially calibrated optical 
data to compute new illumination normalization fac- 
tors until said factors are stabilized. 

41. A method for generating calibrated optical data 
that characterizes a population of cells, the method 

25 comprising: 

(a) scanning a population of calibration particles 
comprising a known constant parameter with a 
beam along a scan line to produce digital cali- 
bration data samples, different said calibration 

30 samples representing different locations along 
said scan line; 

(b) storing said digital calibration data; 

(c) locating a calibration particle within said pop- 
ulation; 

35 (d) defining a neighborhood around the located 
particle, said neighborhood containing optical 
calibration data derived from said stored digital 
calibration data samples within said neighbor- 
hood; 

40 (e) determining the location of a maximum digi- 
tal data calibration sample within said neighbor- 
hood and recentering the neighborhood with re- 
spect to said maximum digital calibration data 
sample; 

45 (f) summing all optical calibration data values in 
a recentered neighborhood to generate an op- 
tical neighborhood calibration value for that 
neighborhood; 

(g) processing said optical calibration data to 
so generate an illumination normalization factor to 

correct optical data for illumination differences 
as a function of particle location along the scan 
line; 

(h) scanning the cell population with a beam 
55 along said scan line to produce digital data 

samples, different said samples representing dif- 
ferent locations along said scan line within the 
cell population; 
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(i) storing said digital data; 

(j) deriving optical data from said stored digital 

data; and 

(k) correcting optical data for each cell at dif- 
ferent locations along the scan line by multiply- 5 
ing the optical data for that cell by the cor- 
responding illumination normalization factor. 

42. The method of claim 41, wherein said process- 
ing of said optical calibration data comprises: 

(a) determining the most frequently occurring io 
optical calibration neighborhood value to define 

a mode optical calibration value; 

(b) determining a set of average optical calibra- 
tion neighborhood values of particles as a func- 
tion of particle location along the scan line; and is 

(c) computing an array of correction coefficients 
equal to the ratio of the mode optical calibration 
value to the average optical calibration value for 
each location along the scan line to generate an 
illumination normalization factor. 20 

43. The method of any one of claims 1, 14, 17, 22, 
35, or 41, wherein the location of each cell found is 
recorded and stored. 

44. The method of any one of claims 1, 14, 17, 22, 

35, or 41, wherein the time corresponding approxi- 25 
mately to the time the digital data was acquired is 
recorded and stored as synchronous digital time 
data, whereby each sample has a corresponding 
digital time point. 

45. The method of any one of claims 1, 14, 17, 22, 30 
or 35, wherein said cell locating step comprises 
comparing 8 said stored digital data to a preselec- 
ted threshold. 
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